Multifractal spectra of mean first-passage time distributions in disordered chains 



Pedro A. PurjQ 

Facultad de Matemdtica, Astronomia y Fisica, 

Universidad National de Cordoba, 
Ciudad Universitaria, 5000 Cordoba, Argentina 

Manuel 0. Cacere^j] 

Centro Atomico Bariloche, Instituto Balseiro, and 
Consejo National de Investigaciones Cientificas y Tecnicas, 
8400 San Carlos de Bariloche, Rio Negro, Argentina 
(Dated: Physical Review E - Received 24 February 2003) 

The multifractal characterization of the distribution over disorder of the mean first-passage time 
in a finite chain is revisited. Both, absorbing-absorbing and reflecting-absorbing boundaries are 
considered. Two models of dichotomic disorder are compared and our analysis clarifies the origin of 
the multifractality. The phenomenon is only present when the diffusion is anomalous. 

PACS numbers: 05.40.Fb, 05.60.Cd, 89.75.Da, 02.50.Ng 



I. INTRODUCTION 

In the past two decades, a great effort has been devoted 
to the study of diffusion and transport in disordered me- 
dia by models based on random walks Q,Q,E|- Basically, 
the dynamics of the system can be described by a master 
equation for a particular probability distribution. Two 
alternative ways are usually employed. The first one is 
based on the probability P(r, t) that the walker is on site 
r at time t when starting from the origin at t = 0. The 
second way consists of analyzing the statistics of the exit 
time from a given region. In both cases, there are ex- 
act enumeration techniques which enable us to calculate 
the corresponding observable in all the possible disorder 
configurations. Thus, we can numerically reckon the dis- 
tribution over disorder of the probability -PjVji), or the 
moments {P q (r,t)) of this distribution 0,1111 IS SHU. 
In an analogous manner, we can numerically evaluate the 
distribution of the mean first-pas sag e time (MFPT) over 
disorder, or its moments |Tot TlTL f]~3L [T^. f]~5L fl6 | ■ 



M(q) 



T q p(T) dT , 



(1) 



where p(T) is the distribution over disorder and q is a 
real number not necessarily integer. The crucial point 
of diffusion in disordered media is that transport can 
be anomalous, i.e., the mean square displacement of a 
random walk scales with time as (r 2 ) = Dt 2 l dw , where 
d w > 2 is the anomalous diffusion exponent. A relevant 
property of the anomalous diffusion is that it leads to 
broad distributions. Their moments cannot be described 
by a single exponent but an infinite hierarchy of expo- 
nents is needed to characterize them H ^jj . This fact 
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enables us to study the distributions over disorder with 
the multifractal formalism 0, 0] . 

Random walks on random fractals (such as the infi- 
nite percolation cluster at criticality) are processes with 
anomalous diffusion, which are characterized by logarith- 
mically broad distribution functions reflecting an under- 
lying multifractal structure jlH HHHITcH] . An impor- 
tant issue in these systems is the nature of the multifrac- 
tality based on the dynamical process or in the steady 
state (as the voltage drop distribution in percolation). 
For the voltage drop distribution, the multifractal be- 
havior is well established for all values of q whereas 
for the dynamical diffusive process the multifractality 
only appears in a range of values of q 8]. Another 
property of these processes is that while the multifrac- 
tal behavior appears in the distribution of the proba- 
bility P(r, t), it is not present in the mean number of 
distinct sites visited by a particle diffusing on the per- 
colating cluster |19| . Another dynamical process with 
"ultra-anomalous" slow motion is the Sinai model [2(j 
for diffusion in a linear chain in the presence of random 
fields. For this model, the mean square displacement of 
a random walk follows asymptotically (r 2 ) w ln 4 i. Sev- 
eral studies H IH El El El have established the 
multifractal properties of the model. However, the ori- 
gin of the multifractality in this process is yet an open 
question. 

In this work, we revisit the Sinai model and we address 
the last question. Particularly, we compare the behavior 
of the MFPT distributions over disorder and its moments 
for two processes with dichotomic disorder. One of them 
is the Sinai model and the other is a nonanomalous biased 
random walk in a finite disordered chain. The outline of 
the paper is as follows. In Sec. [HI we present expressions 
for the MFPT for a given realization of the (quenched) 
disorder in the chain. The description of the models of 
disorder, employed in our studies, is given in Sec. lIIII The 
moments of the MFPT distribution over disorder diverge 
with the system's size. These divergences are character- 
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ized by the scaling exponents in S ec. I IV1 where we 
employ the multifractal formalism 0, Il8| to calculate 
the exponents r(q) of the corresponding partition func- 
tion. Also, the generalized Renyi dimensions D(q), and 
the spectra f(a) are given in that section. Finally, in 
Sec. we briefly summarize the results of our work. 



II. MEAN FIRST-PASSAGE TIME IN 
QUENCHED DISORDERED CHAINS 

We consider the continuous time dynamics of a random 
walk on a discrete one-dimensional lattice with nearest 
neighbor hopping. The walker jumps from site n to site 
n + 1 with transition probability per unit time or 
to site n — 1 with transition probability per unit time 
w~ . We are concerned with the exit time of the walker 
from the finite interval D = [—ill, L] on the chain, with 
at least one absorbing end. The average of the survival 
time until the absorption, over realizations of the random 
walk, is the MFPT. Recently, we have obtained a general 
exact expression for the MFPT for a fixed set of transi- 
tion probabilities {wf} 22j. Let T n denote the MFPT if 



the walker initially began at site n E D. For an interval 
with both absorbing extremes (w~^, M+l -. = w~ L+1 = 0) 
we get 



e n 



T„ = 



k=-M j=-M 3 

L k - 

i+ £ n % 

k=-M j = -M 3 

L L-l 



X 



5w W k k^M W k W t , 

n—1 _^ n — 2 ^ n—1 i ^ — 

e z+ + 1! 7 e n 3 

k fc=_M W k k=-M W k i=fc+lj=fc+l W J 



(2) 



In this work we fix the starting point n — (defining 
T = Tq) and we consider two possible boundary condi- 
tions. The first one is the interval [—L,L] with absorbing- 
absorbing (AA) extremes. The total number of sites, N s , 
in the interval is 2L + 1. The MFPT for this case is given 
by Eq. (J2Jl taking M = L. The second case is the in- 
terval [0, L] with reflecting-absorbing (RA) ends. Here, 
N s = L + 1 and from Eq. ©, taking n = 0, M = 0, and 



0, we immediately obtain for the MFPT [2J 



L-l 



T = 



Ei + Ei e n 5 



fc=0 W & i=k+lj=k+l W 3 



(3) 



The effect of the reflecting boundary is the striking sim- 
plification of the structure of the equation for the MFPT. 



This fact leads us to consider also the AA boundary con- 
ditions in the problem of multifractality of the MFPT 
distribution over disorder. 

In finite discrete systems, we can enumerate all the 
configurations N of disorder. We denote by T/f, the 
MFPT for the zth realization of the quenched disor- 
der. We stress that all the configurations of disor- 
der are equally probable. However, the resulting val- 
ues T {l) are distributed by p{T), the MFPT distribution 
over disorder. Thus, we can compute exactly the mo- 
ments of the MFPT distribution, given the set of values 
{T(i\, i = 1, . . . ,jV}, from the definition given by Eq. lfT)l. 



M(q) 



1 



AT 

En 



(4) 



In particular, the MFPT averaged over disorder results 
in M(l). For dichotomic models of disorder, in a chain 
with N a sites, there are j\f = 2 Ns possible realizations 
of the random lattice, which can be easily enumerated. 
The set of values {T(i), i = 1, ...,Af} can be exactly 
calculated employing expression @ for the AA boundary 
conditions, or by Eq. © for the RA extremes. 



III. MODELS OF DISORDERED CHAINS 

Now, we assume that the hopping probabilities are 
strictly positive random variables, chosen independently 
from site to site and identically distributed. Additionaly, 
we admit that the site transition probabilities are not 
necessarily symmetric in the sense that w1~ ^ wj . Thus, 
we can incorporate the effects of bias in to the chain by 
external fields. 

We select the first distribution for the transition prob- 
abilities in such a way that they satisfy the Sinai condi- 
tion [2(| , namely, the random variable In (w'j' / wj ) has 
zero mean and finite variance a 2 . Thus, we consider a 



dichotomic model by defining w- = 1 



and pre- 



scribing wj~ to be equal to 1/2 ± e with equal probabili- 
ties. The parameter e measures the strength of disorder, 
and can take values between and 1/2. e = corre- 
sponds to a simple homogeneous random walk, and for 
e > we get a disordered random walk with local bias 
(random field). It is easily verified that the above pre- 
scriptions satisfy the Sinai condition. Particularly, we 
obtain a 2 = ln 2 7(e), where 7(e) = (1 + 2e)/(l - 2e). 
In the asymptotic limit e — » 1/2, the variance diverges. 
The MFPT averaged over disorder for the dichotomic 
Sinai model diverges as B(e) " for N s — ► 00, where 
(3(e) = (1 + 4e)/(l — 4e) [23, whereas the typical value 
of the first-passage time, defined as exp((lnT)), diverges 

slower, as exp (^t^itN s /2^ Here, the brackets (• • • ) 

denote the average over disorder. Therefore, the distribu- 
tion of the MFPT over the Sinai disorder has a power-law 
tail [H. 
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FIG. 1: Numerical evaluation of Tmax as a function of N 3 
for chains with the AA ends and L from 1 to 13. 



It is easily seen that the largest value of the MFPT, 
Tmax, is obtained for the RA ends when at all the sites, 
the right jump transition is 1/2 — e, and the left jump 
transition is 1/2 + e. In the asymptotic limit of N s — > oo, 
we obtain Tmax ~ T^ 8 G3- For the AA extremes, T ma x 
cannot be easily evaluated. Thus, instead of an analytical 
expression, in Fig. ^ we show a numerical calculation of 
Tmax for several values of e. We see that asymptotically 
log(Tmax) ~ N s , as in the RA case. 

In Fig. |3 we show the histogram of the distribution 
p(T) over the dichotomic Sinai disorder. The figure was 
constructed computing Eq. J5J for M = L = 7 (N s = 
15), and using all the possible configurations of disorder. 
A related figure, for the histogram of p(T) for chains 
with the RA extremes, can be seen in Ref. 1141 . The 
distribution p(T) has previously been studied 25] for the 
RA boundary conditions and results broad. 

Our second random biased model for the transition 
probabilities is defined by w~ = — e and prescribing 



to be equal to 1/2 or 3/2 with equal probabilities. In 
this case, the parameter e measures the strength of the 
bias, and can take values between and 1/2. In the limit 
of e = 0, we get a disordered symmetrical random walk. 
This dichotomic model corresponds to a class of weak dis- 
order with global bias. The quantities /?& = ^(l/w+) fe ^ 

result finite for all k > 1. Therefore, the model does 
not present anomalous diffusion For the MFPT av- 
eraged over disorder with the A A boundary conditions, 
we obtain up to first order in e [2^, [2(j , 



M(l) 



(L + l) 2 
Wi 1 



(5) 



FIG. 2: Histogram of p(T) for a chain with L = 7, the AA 
boundary conditions, and the Sinai disorder with e = 0.25. 
We are using the rescaled variable In T/ In Tmax and uniform 
buckets of width equal to 2 9 . Tmax = 9824. 



The asymmetry in the hopping transitions links the 
strength of the bias with the fluctuation of the disorder, 
defined by T = (ft — 0i) / f3f- In our particular case, we 
get Pi — 4/3, and T — 1/4. On the other hand, for the 
RA extremes, we obtain up to first order in e |22|. 



M(l) 



(L + 1)(L + 2) 
2/3^ 



(6) 



Strikingly, for these boundary conditions, the fluctuation 
of the disorder is not present in the averaged MFPT. This 
fact relies on the difference in the structures of Eqs. 
and J3J)- In the limit of L — * oo, M(l) diverges as L 2 
independently of the boundary conditions. 

In Fig. [31 we show the histogram of the corresponding 
distribution p(T) over dichotomic weak disorder. Here, 
we do not obtain a broad but a more localized distribu- 
tion. We must stress the difference in the scales used to 
construct the plots in Figs. |21 and |3J This fact is the first 
indication of the different nature in the distribution over 
disorder between both models. Both models were con- 
structed to get = 1, and the main difference among 
them is in the role of the parameter e. This parameter 
controls the disorder in the Sinai model and regulates the 
bias in the second model. Thus, in the Sinai model the 
bias is local, i.e., the direction of the bias is randomly 
drawn in each site, whereas in the second model, we are 
considering a global bias field which points to the right in 
all sites. Assuming that the jump probabilities per unit 
time involve the Arrhenius factor (see Ref. 26]), the land- 
scape for the particle potential is schematically sketched 
in Fig. 2| for both models. 
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FIG. 3: Histogram of p(T) for a chain with L = 7, the AA 
boundary conditions, and weak disorder with e = 0.25. We 




are using uniform intervals of width equal to 2 
31.75. 



Tmax 



FIG. 4: Sketches of the particle potential for both kinds of 
disorder and small e: (a) The Sinai model. The difference 
between two consecutive peaks is 2e, with random sign, (b) 
Weak disorder. The difference between one peak and the 
following is always e. 



IV. MULTIFRACTAL ANALYSIS 

In order to characterize the divergences of the moments 
of the MFPT distribution over disorder, we postulate 
that in the limit N s — > oo, the qth moment obeys the 
scaling relation M(q) ~ Af^^ . In order to verify this 
ansatz, we have numerically reckoned log[M(g)] vs N a . 
For chains with the AA ends, we take L from 1 to 13 and 
use Eq. for evaluating the set {Tu\, i = 1, . . . , A/"}, for 
a given L. For chains with the RA boundary conditions, 
we take L from 1 to 20 and use Eq. @ for calculating 
the values To* , for a fixed size of the chain. For all values 
of e and q, and for both models of disorder, with the AA 
or RA extremes, we obtain straight lines. 

The slope of the linear fitting for the Sinai model with 
the AA (RA) ends is plotted in Fig.[51(Fig.[BJ| as functions 
of q, for three values of e. For all e, we obtain that £(q) is 
a monotonically increasing function with £(0) = 0, which 
is due to the normalization condition of the distribution. 
We observe that for large values of e, is a nonlinear 
function; whereas in the limit of weak disorder, be- 
comes a linear function. For the model of weak disorder, 
we obtain = 9q, where 9 > is a decreasing func- 
tion of e. Therefore, for weak disorder, we get a single 
gap exponent for describing the moments of the MFPT 
distribution. The analysis with both kinds of boundary 
conditions leads us to qualitatively similar behaviors for 
the exponents £(g), for both models of disorder described 
in Sec. ED 

Now, let us define the partition function as follows : 
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FIG. 5: Plot of the scaling exponents for chains with 
the AA boundary conditions and the Sinai disorder. 





FIG. 6: Plot of the scaling exponents £(5) for chains with 
the RA boundary conditions and Sinai disorder. 

Again, we postulate that in the limit N s — > 00, the par- 
tition function obeys a scaling relation. Thus, we write 
Z(q) ~ Af~ T ( q \ It is well known that for nonmultifrac- 
tal distributions, the exponents r(q) are linear functions 
on q, namely, r(g) = q — 1. The multifractality appears 
with a nontrivial dependence of the scaling exponents on 
q. For all values of e and q, and for both models of dis- 
order, we obtain straight lines when we plot log[Z(g)] vs 
N s . Therefore, the scaling ansatz is verified. Here, we 
have also used L from 1 to 13 for the AA ends, and L 
from 1 to 20 for the RA extremes. 

In Fig. (Fig. [HJl, we plot r(q) for chains with the di- 
chotomous Sinai model, the AA (RA) ends, and for three 
representative values of e. The functions r(q) satisfy 
t(0) = —1. Additionally, we get r(l) = 0, which is due 
to the normalization condition of Z(q). In the limit of 
e — > 0, r(q) ~ q — 1, and for strong disorder, r{q) exhibits 
a nonlinear dependence. Thus, the multifractality for the 
Sinai model is unambiguously established, independently 
of the boundary conditions. For the model of weak dis- 
order, we obtain r{q) = q — 1 for all values of e, within 
the accuracy of the numerical evaluations. Therefore, we 
have not found multifractality in our model of weak dis- 
order, neither with the AA nor with the RA boundary 
conditions. This result reinforces the idea that the mul- 
tifractal phenomenon is only related to anomalous diffu- 
sion. Both models are based on the same dichotomic rule 
for assigning values to the transition probabilities. More- 
over, in both cases we found (u^) = 1. In spite of these 
similar characteristics, the selected models have quite dif- 
ferent behaviors. The random field resulting from the 
Sinai condition is at the basis of the anomalous diffusion 
and the origin of the multifractal behavior of the MFPT 



FIG. 7: Plot of the scaling exponents r(q) for chains with 
the AA boundary conditions and the Sinai disorder, using L 
from 5 to 13. The solid line corresponds to the linear relation 
q-1. 



distribution over disorder. 

The simple plots obtained in the graphs for the func- 
tions r(q) suggest us that the multifractality is due to 
the presence of a binomial multiplicative process. For 
this process, the mass exponents are 

r (5 )=-^±M), (8) 

where p\ + P2 = 1 ■ This expression immediately satisfies 
the conditions r(0) = — 1 and r(l) = 0, and in the limit 
Pi = Pi = 1/2, results in r(q) = q — 1. For p\ < 1/2, we 
get the limit expressions 

r(q) -> -q for q -> -00 , (9) 

r(q) -> -q ^y for q -> 00 . (10) 

In our case, the parameters p\ and P2 are functions of 
the strength of disorder (e). In Table U the values of the 
fitting of r(g), using the two parameters p\ and pi-, are 
displayed for both kinds of boundary conditions. We only 
obtain a good quality of fitting for small values of e. For 
strong disorder, the condition p\ + P2 = 1 is relaxed and 
we can only fit both asymptotic regimes of the curves for 
q — » ±00. This fact is a strong indication that the nature 
of the multifractal phenomenon is more complex than a 
multiplicative rule. 

A usual characterization of multifractality is the gen- 
eralized Renyi dimension spectrum |18j defined by |27| 



FIG. 8: Plot of the scaling exponents r(q) for chains with 
the RA boundary conditions and the Sinai disorder, using L 
from 6 to 14. The solid line corresponds to the linear relation 
9-1. 



TABLE I: Fitted values for the parameters pi and P2, using 
Eq. ©, for the plots r(q) in Figs. and [HI 



e 


Pi 


AA 

pi 


Pi 


RA 

pi 


0.00 


0.500 


0.500 


0.500 


0.500 


0.10 


0.461 


0.543 


0.425 


0.602 


0.25 


0.395 


0.700 


0.295 


0.878 


0.35 


0.336 


0.825 


0.173 


0.964 



m - T[q) 



9-1 



(11) 



In Fig. [5] (Fig. 1 1 (J|) we depict the generalized dimension 
spectra for the same conditions and values of e as given 
by Fig. [7| (Fig. [8J. We find that D(q) is a monotonically 
decreasing function of q, and satisfies D(0) — 1, which is 
the dimension of the support of the distribution. 

Taking the Legendre transform of r(q), we can obtain 
the multifractal spectrum /(a) = —t(q) + qa, where 
the Lipschitz-Holder exponent a is the derivative of r(q) 
with respect to q. Figure ITD (Fig. IT^|) shows the f(a) 
spectra for the same conditions and values of e as given 
by Fig. (Fig. 0. For strong disorder, the curve f(a) 
becomes broad, whereas for e — > the spectrum collapses 
to the point (1,1). As is known, the maximum value of 
/ is the fractal dimension of the support of the measure. 

In this work, we have chosen the scaling parameter A/" 
to characterize the system size. In Ref. 14], a similar 



FIG. 9: Generalized dimensions D(q) for chains with the AA 
boundary conditions and the Sinai disorder. 
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FIG. 10: Generalized dimensions D(q) for chains with the 
RA boundary conditions and the Sinai disorder. 



multifractal analysis for the Sinai disorder and the RA 
boundary conditions was based on the scaling parameter 
Tmax- We have seen in Sec. IIIII that Tmax oc A/". How- 
ever, this choice of the scaling parameter leads to spuri- 
ous results, such as the variation of D(0) with e. This 
phenomenon is an artifact of the mathematical selection 
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of the scaling parameter, and not a characteristic prop- 
erty of the system under analysis. On the other hand, for 
one-dimensional systems, we expect that the dimension 
of the support of the multifractal measure is equal to 1 . 

Finally, we derive the relation between the exponents 
£(</) and r(q). From the definitions given by Eqs. Q 
and J7J), we inmediatly obtain 



Z(q) = M(q) 
w M (l)o ' ' 



(12) 



Using the scaling ansatz in the last equation and taking 
logarithms, we get 



r(q) 



(13) 



Particularly, the last relation satisfies inmediately that 
t(0) = — 1 and r(l) = 0. From the result quoted in 
Sec. IIIII for chains with the RA ends, we found £(1) = 
log 2 /?(e). 



a 



FIG. 11: Multifractal spectra f(a) for chains with the AA 
boundary conditions and the Sinai disorder. 
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FIG. 12: Multifractal spectra f(a) for chains with the RA 
boundary conditions and the Sinai disorder. 



V. CONCLUDING REMARKS 

In this work, we have considered the distribution of the 
MFPT over two classes of disorder, the Sinai and another 
dichotomic model with global bias and nonanomalous dif- 
fusion. Our results confirm us that the multifractality is 
related only to anomalous diffusion. 

The multifractal behavior in the MFPT distribution 
over the Sinai disorder is not a consequence of the mul- 
tiplicative structure of Eq. (J2J or © (the expressions for 
the MFPT for a fixed realization of disorder), as was 
stated in a previous report 14]. The multifractality is 
an inherent attribute of the strong disorder in the Sinai 
model, and is well established for both kinds of boundary 
conditions. Moreover, the multifractal signature found 
in the spectra obtained suggests us that the origin of the 
phenomenon is more complex than a binomial multiplica- 
tive process. 



Acknowledgments 

This work was partially supported by grant from "Se- 
cretaria de Ciencia y Tecnologia de la Universidad Natio- 
nal de Cordoba" (Code: 05/B160, Res. SeCyT 194/00). 
P.A.P. thanks Centro Atomico Bariloche for its hospital- 
ity during his stay. P.A.P. is also grateful to Peter Pfcifer 
for a useful discussion. 



[1] S. Alexander, J. Bernasconi, W. R. Schneider, and R. Or- 

bach, Rev. Mod. Phys. 53, 175 (1981). 
[2] J. W. Haus and K. W. Kehr, Phys. Rep. 150, 263 (1987); 

S. Havlin and D. Ben-Avraham, Adv. Phys. 36, 695 



(1987); J. P. Bouchaud and A. Georges, Phys. Rep. 195, 
127 (1990). 

[3] S. Havlin and A. Bunde, in Fractals and Disordered Sys- 
tems, edited by A. Bunde and S. Havlin (Springer- Verlag, 



8 



Berlin, 1991), Chap. 3. 
[4] A. Bunde, H. Harder, S. Havlin, and H. E. Roman, J. 

Phys. A 20, L865 (1987). 
[5] H. E. Roman, A. Bunde, and S. Havlin, Phys. Rev. A 

38, 2185 (1988). 
[6] A. Bunde, S. Havlin, and H. E. Roman, Phys. Rev. A 

42, 6274 (1990). 
[7] H. E. Roman, Physica A 191, 379 (1992). 
[8] E. Eisenberg, A. Bunde, S. Havlin, and H. E. Roman, 

Phys. Rev. E 47, 2333 (1993). 
[9] E. Eisenberg, S. Havlin, and G. H. Weiss, Phys. Rev. 

Lett. 72, 2827 (1994). 
[10] N. Masic and Z. V. Djordjevic, Physica A 167, 560 

(1990). 

[11] C. Van den Broeck, J. Stat. Phys. 65, 971 (1991). 

[12] K. P. N. Murthy, S. Rajasekar, and K. W. Kehr, J. Phys. 
A 27, L107 (1994). 

[13] T. Wichmann, A. Giacometti, and K. P. N. Murthy, 
Phys. Rev. E 52, 481 (1995). 

[14] K. P. N. Murthy, K. W. Kehr, and A. Giacometti, Phys. 
Rev. E 53, 444 (1996). 

[15] K. P. N. Murthy, A. Giacometti, and K. W. Kehr, Phys- 
ica A 224, 232 (1996). 

[16] K. Kim, J. S. Choi, and Y. S. Kong, J. Phys. Soc. Jpn. 
67, 1583 (1998); K. Kim, G. H. Kim, and Y. S. Kong, 
Fractals 8, 181 (2000). 



[17] T. C. Halsey, M. H. Jensen, L. P. Kadanoff, I. Procaccia, 
and B. I. Shraiman, Phys. Rev. A 33, 1141 (1986). 

[18] J. Feder, Fractals (Plenum Press, New York, 1988). 

[19] L. K. Gallos, P. Argyrakis, and K. W. Kehr, Phys. Rev. 
E 52, 1520 (1995); K. P. N. Murthy, L. K. Gallos, P. Ar- 
gyrakis, and K. W. Kehr, ibid. 54, 6922 (1996). 

[20] Y. G. Sinai, Theor. Probab. Appl. 27, 256 (1982). 

[21] R. Harisha and K. Murthy, Physica A 287, 161 (2000). 

[22] P. A. Pury and M. O. Caceres, J. Phys. A 36, 2695 
(2003). 

[23] K. P. N. Murthy and K. W. Kehr, Phys. Rev. A 40, 2082 

(1989) ; 41, 1160(E) (1990). 

[24] S. H. Noskowicz and I. Goldhirsch, Phys. Rev. Lett. 61, 
500 (1988). 

[25] K. W. Kehr and K. P. N. Murthy, Phys. Rev. A 41, 5728 

(1990) ; S. H. Noskowicz and I. Goldhirsch, ibid. 42, 2047 
(1990). 

[26] P. A. Pury and M. O. Caceres, Phys. Rev. E 66, 021112 
(2002). 

[27] Note that our definition of r(q) is equal to minus the mass 
exponents used in Ref. [l8j|. Usually, the mass exponents 
are defined in the limit of a scaling parameter going to 
zero. In our case, the scaling relations hold for the param- 
eter M — > oo. Therefore, we use the denominator (q — 1) 
in the definition of D(q). 



